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ABSTRACT 

The origin of massive field stars in the Large Magellanic Cloud (LMC) has long been an enigma. The recent measurements of large 
offsets (~ 100 kms~') between the heliocentric radial velocities of some very massive (02- type) field stars and the systemic LMC 
velocity provides a possible explanation of this enigma and suggests that the field stars are runaway stars ejected from their birthplaces 
at the very beginning of their parent cluster's dynamical evolution. A straightforward way to prove this explanation is to measure the 
proper motions of the field stars and to show that they are moving away from one of the nearby star clusters or OB associations. This 
approach is, however, complicated by the long distance to the LMC, which makes accurate proper motion measurements difficult. 
We used an alternative approach for solving the problem (first applied for Galactic field stars), based on the search for bow shocks 
produced by runaway stars. The geometry of detected bow shocks would allow us to infer the direction of stellar motion, thereby 
determining their possible parent clusters. In this paper we present the results of a search for bow shocks around six massive field stars 
that have been proposed as candidate runaway stars. Using archival Spitzer Space Telescope data, we found a bow shock associated 
with one of our programme stars, the 02 V((f*)) star BI237, which is the first-ever detection of bow shocks in the LMC. Orientation 
of the bow shock suggests that BI 237 was ejected from the OB association LH 82 (located at = 120 pc in projection from the star). A 
by-product of our search is the detection of bow shocks generated by four OB stars in the field of the LMC and an arc-like structure 
attached to the candidate luminous blue variable R81 (HD 269128). The geometry of two of these bow shocks is consistent with the 
possibility that their associated stars were ejected from the 30 Doradus star-forming complex. We discuss implications of our findings 
for the problem of the origin of runaway stars and the early dynamical evolution of star clusters. 

Key words. Stars: kinematics and dynamics - stars: individual: BI237 - stars: individual: HD 269128 - open clusters and associa- 
tions: individual: LH 82 - open clusters and associations: individual: R136 (HD 38268) 



o 

6 

CO 

<N 

> : 

' 1. Introduction low significance (Gvaramadze & Bomans [2008bt Gvaramadze 

, et al. I2Q 1 Qbl > . It is therefore tempting to search for bow shocks 

O ' Although the majority of massive stars are situated in their par- around field OB stars in me Lar S e Magellanic Cloud (LMC) 

■ ent clusters and OB associations, a significant population of where accurate P ro P er motlon measurements are difficult, while 

O i young massive stars exists in the field, some of which are sep- bow shocks can stl11 be resolved wlth modern lnfrared tele " 

<0 arated by hundreds of parsecs from known clusters and OB sc °pes. 

, associations (Garmany, Conti & ChiosiEMI Garmany 1990} In thls P a P er we P resent the results of a search for bow 

^ ■ Massey & ContiUMJ Massey et al.QSS). Some Galactic field shocks ln the LMC usln 8 S P ltzer S P ace Telesc °P e data - 

.£ : stars have high measured peculiar (either radial or transverse) ° ur P nme S oal was to detect bow shocks Produced by iso- 

X ■ velocities and are therefore most likely runaway stars ejected lated ' ver y masslve stars that have previously been qualified 

H ' from a cluster (Blaauw 1961 ; Gies & BoltonM StoneHSSJ as runawa y s on the basls of their lar S e P ecullar radlal velocl " 

. P. , Zinnecker 2003 ). A straightforward way to prove the runaway tles < Sect 2 >' We ^covered a bow shock associated with one 

nature of the field OB stars is to use the available kinematic data of these stars ' the 02 V « f *» star BI 231 ' A by-product of our 

on these stars to back-trace their orbits to parent clusters (e.g. search 18 Section of bow shocks produced by several other iso- 

Hoogerwerf, de Bruijne & de Zeeuw 2001 ). Schilbach & Roser lated 0B stars ' and two of these stars m located around the 

(2008 ) make extensive use of this approach to show that most 30 Doradus nebula < SecL 3 >' Implications of our findings for the 

Galactic field OB stars are formed in clusters. Alternatively, the P roblem of the on g ln of mnawa y stars and the earlv d y naml " 

runaway nature of the field OB stars can be proved via detec- cal evolution of star clusters are discussedm Sect. 4. We use a 

tion of their bow shocks - the natural attributes of supersoni- dlstance °{ 50 k P c for the LMC (Gibson[200S so that 1 corre- 

cally moving objects (e.g. Baranov, Krasnobaev & Kulikovskii s P on ds to - 14 pc. 
[19711 Van Buren & McCray IT9881 . The geometry of detected 

bow shocks would allow one to infer the direction of stellar mo- ^ Verv massive field Stars as runaways 

tion (Van Buren, Noriega-Crespo & Dgani 1995), thereby deter- 
mining the possible parent clusters even for those field OB stars The study of the massive star population in the LMC by Massey 
whose proper motions are still not available or measured with a et al. (1995 ) has shown that a large number of young very mas- 



1 



V.V. Gvaramadze et al.: Massive runaway stars in the Large Magellanic Cloud 



sive (02-type) stars is located at ~ 100 - 200 pc in projection 
from star clusters and OB associations. This finding was inter- 
preted as indicating that the field can produce stars as massive as 
those born in clusters (Massey et al. 11995b . An obvious alterna- 
tive to this interpretation is that the massive field stars were actu- 
ally formed in a clustered environment and subsequently ejected 
from their birth sites via dynamical processes (e.g. Clarke & 
Pringle [T992l Pflamm-Altenburg & Kjoupa [2006l ); i.e., the field 
massive stars are runaway stars (Walborn et al. 2002 ; Brandl et 
al. 120071 1. The large separations from the possible parent clus- 
ters and the young (~ 2 Myr) ages of the very massive field 
stars imply that their (transverse) velocities should be as high as 
~ 50 - 100 km s 1 (Walborn et al. [20021 . The large offsets from 
the parent clusters and the high peculiar velocities are not un- 
usual for Galactic massive runaway stars. For example, the 04If 
star BD+43 3654 ejected from the CygOB2 association is lo- 
cated at about 80 pc from the core of the association (Comeron & 
Pasquali f2007t Gvaramadze & Bomans 2008a), while the helio- 
centric radial velocity of the star is offset by 56 km s _1 from the 
systemic velocity of Cyg OB2 (Kobulnicky, Gilbert & Kiminki 
2010), which with the transverse peculiar velocity of the star of 
^ 40 kms~' (Comeron & Pasquali 2007 ) corresponds to the to- 
tal peculiar velocity of ^ 70 km s _1 . 

The runaway interpretation of the massive field stars in the 
LMC received strong support after the discovery that some of 
them have high peculiar radial velocities, much greater than the 
systemic velocity of the LMC (Massey et al. 120051 Evans et al. 
125061 120T01 Earlier, Nota et al. ( fT99lj ) and Danforth & Chu 
(120011 ) had found that the systemic velocity of the candidate lu- 
minous blue variable SI 19 (HD 269687) and its circumstellar 
nebula is ~ 100 km s lower than that of the LMC, so they sug- 
gested that SI 19 could be a runaway star. 

To prove the runaway nature of the field massive stars in 
the LMC unambiguously, it is necessary to ascertain their par- 
ent clusters. The proper motion measurements cannot help solve 
this problem in the near future. At the distance of the LMC, 
the transverse velocity of ~ 50 - 100 kms -1 corresponds to a 
proper motion of ~ 0.2 - 0.4masyr _I , which is too low to be 
measured with high confidence using the ground-based obser- 
vations. The simpler solution is to infer the direction of stellar 
motion via the geometry of bow shocks produced by runaway 
stars (Gvaramadze & Bomans 2008b ) or, in the case of evolved 
massive runaway stars, through the asymmetry in the bright- 
ness distribution of associated circumstellar nebulae (Danforth 
& Chu 2001 1 Gvaramadze et al. 2009). It is worth noting that two 
very massive Galactic runaway stars, BD+43° 3654 and ACep 
(06I(n)f; Walborn 1973), are both associated with spectacular 
bow shocks (Comeron & Pasquali 2007; Van Buren et al. U995b . 
One can therefore expect that some of the field stars in the LMC 
will manifest themselves in such secondary attributes of high- 
velocity runaway stars. 

3. Search for bow shocks in the LMC 

To search for bow shocks, we selected four isolated massive stars 
with high peculiar radial velocities (Massey et al. 2005 ; Evans et 
al. [2006 2010) and added two isolated 02-type stars, BI253 
(02 V((f* ))) and Sk-68°137 (02III(f*)), which Walborn et 
al. ((2002) suggests are runaways owing to their large separa- 
tion from their plausible birthplace in the central cluster, R136 
(HD 38268), of the 30Doradus nebula. The details of these stars 
(listed in order of their RA) are summarized in TableQ] For the 
first four stars, we give their peculiar radial velocities, while the 
transverse velocities are listed for the remaining two, inferred 



Table 1. Summary of candidate runaway stars in the LMC. 



Star 


Spectral type 


v ( km s 1 ) 


N 11-026 


02.5 III(f*) 


=s 35 {a> 


Sk-67°22 


02 If* 


a 150<» 


BI237 


02V((f*)) 


~ 120m 


30 Dor 016 


02 Ill-If* 


=! 85 (c) 


BI253 


02V((f*)) 


~ 55» 


Sk-68°137 


02III(f*) 


~ 100 (</) 



Notes. (fl) Evans etal.l2006l (b) Massey et al. l2005l (c) Evans et al. l20T0l 
<rf) Walborn etal. l2002l 



under the assumption that both stars were ejected ~ 2 Myr ago 
from R136 (Walborn et al.|2002J. 

From our experience in the search for bow shocks produced 
by OB stars ejected from Galactic star clusters (Gvaramadze & 
Bomans 2008a b; Gvaramadze et al. 1201 Obi ) using the archival 
data from the Midcourse Space Experiment (MSX) satellite 
(Price et al. 120011 ) and the Spitzer Space Telescope (Werner et 
al. 1200 4). we know that the bow shocks are visible mostly in 
21.3 pm (MSX band E) images and 24 pm images obtained with 
the Multiband Imaging Photometer for Spitzer (MIPS; Rieke et 
al. 120041 ). The resolution of Spitzer 24 pm images (~ 6") is three 
times better than those of the MSX, so that in the search for bow 
shocks in the LMC we utilized the MIPS data alone. 

The typical (transverse) size of bow shocks generated by 
Galactic OB stars (i.e. the extent of a bow shock in the direction 
perpendicular to the vector of the stellar motion) is several par- 
sees; e.g. the size of the bow shocks associated with the above- 
mentioned massive runaway stars, BD+43 3654 and ACep, is 
5.0 and 2.3 pc, respectively. If placed at the distance of the 
LMC, these bow shocks will have an angular size of 10"-20", 
which is comparable to or several times greater than the angular 
resolution of the MIPS 24pm images. Thus, the bow shocks in 
the LMC can be resolved with the Spitzer imaging data! 

Visual inspection of MIPS 24 pm imagesj of fields contain- 
ing our programme stars revealed a bow shock associated with 
only one of them, namely BI237. The non-detection of bow 
shocks around the remaining five programme stars is consis- 
tent with the observational fact that only a small fraction (^ 20 
per cent) of runaway OB stars produce (observable) bow shocks 
(Gvaramadze & Bomans 2008b and references therein). 

FigureQ]gives an overview of the region northwest of BI 237 
with two associations, LH 88 and LH 82, whose centres are sep- 
arated in projection by 65 and 120 pc from the star. (The ap- 
proximate boundaries of the associations are indicated by dashed 
circles; Bica et al. 1999 ) The orientation of the bow shock 
generated by BI237 (see Fig.[2]) suggests that the more likely 
parent association of the star is LH 82. LH 82 contains another 
very massive star, Sk -67° 21 1 (02 III(f*); Walborn et al. 120051 . 
Assuming that BI 237 was indeed ejected from LH 82 and given 
the young (~ 2 Myr) age of the star, one finds that its transverse 
velocity should be ~ 60 km s _1 if the star escaped from the core 
of LH 82 soon after the birth or higher if the ejection event oc- 
curred later on, so that the total peculiar velocity of the star is 
> 130 km s _1 . The angular size of the bow shock of 20" corre- 
sponds to the linear size of 4.8 pc, i.e., a figure typical of mas- 
sive runaway stars (see above). Using these estimates, one can 
constrain the number density of the ambient interstellar medium, 

1 The images, obtained in the framework of the Spitzer Survey of the 
Large Magellanic Cloud (Meixner et al. 2006 ), were retrieved from the 
NASA/IPAC Infrared Science Archive (http://irsa.ipac.caltech.edu). 
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Table 2. Details of bow shock-producing stars. 



N 



LH 82 



J 



LH 88 







Sk -67 211 



o 



Bl 237 



2 arcmin 



Fig. 1. MIPS 24 yum image of the associations LH 82 and LH 88 
(indicated by dashed circles). The position of the 02V((f*)) star 
BI237 and its bow shock are marked by a black solid circle, 
while the position of the 02III(f*) star Sk-67°211 in LH82 
is indicated by a white circle. At the distance of the LMC, 1' 
corresponds to 14 pc. 



30" 



Fig. 2. Left: MIPS 24 fim image of the bow shock associated with 
the 02V((f*)) star BI237. The position of BI237 is marked by 
a circle. Right: 2MASS J band image of the same field. 



«ism- For the characteristic (transverse) size of a (parabolic) bow 
shock of ~ 2 y/2Ro, whereto = (Mv<x,/4-^pisMvi) 1 ^ 2 is the stand- 
off distance of the bow shock, M and Voo are the mass-loss rate 
and terminal velocity of the stellar wind, pism = 1 -4»Jh«ism is 
the density of the interstellar medium, ran the mass of the hy- 
drogen atom, and v„ (> 130 km s _1 ) the velocity of the star rel- 
ative to the ambient medium, and using the wind parameters of 
BI23 7,M = 7.8xl(T 7 Moyr 1 and = 3400 kms" 1 (Mokiem 
et al. 120071 1. one finds «ism < 0.1 cm 3 , i.e., a reasonable num- 
ber. 

A by-product of our search is the detection of bow shocks as- 
sociated with four other OB stars in the field of the LMC (Fig[3j). 
The details of these stars are given in Table[2] The spectral types 
of the stars were found using the VizieR catalogue access tooQ. 
The last column gives the possible birthplaces of the stars. 

The bow shock produced by Sk -69° 206 is located at 17' 
(^ 240 pc) to the west of R136 - the central cluster of the 
30Doradus nebula (Fig.|4|. The orientation of the bow shock 
is consistent with the possibility that the associated star was 
ejected from the 30Doradus nebula. Assuming that the ejection 
event occurred ~ 2 Myr ago, one finds the transverse velocity 
of Sk-69° 206 of 120 kms -1 , which is comparable to the pe- 



Star 




Spectral type 


Association 


Sk- 


66° 16 


09.71b"" 


KMHK268 


Sk- 


68° 86 


Bl:'*' 


[SL63] 495 


Sk- 


69° 206 


B2: w 


30 Doradus 


Sk- 


69° 288 


B0.5 (W 


30Doradus ? 



Notes. <fl) Evans et al. [20061 <*> Rousseau et al. [T978l 
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http://webviz.u-strasbg.fr/viz-bin/VizieR 



Fig. 3. MIPS 24 yum images of bow shocks associated with (a) 
Sk-69° 206, (b) Sk-69° 288, (c) Sk-66° 16, and (d) Sk-68° 
86. The orientation and the scale of the images are the same. 



culiar velocities of the candidate runaway stars in the LMC (see 
Table©. 

The bow shock produced by Sk -69° 288 is situated (at least 
in projection) within the association LH 113. The geometry of 
the bow shock suggests that Sk -69° 288 is moving away from 
the 30 Doradus nebula, which is located at ^ 21' (^ 300 pc) 
to the west of the star (Fig.|4j». Although we have no arguments 
against the association between Sk-69° 288 and LH 113, one 
cannot exclude the possibility that the actual birthplace of the 
star is the 30 Doradus nebula. In this connection, it is worth 
noting that the 05 V star ALS 19631 (Hanson [2003) was sug- 
gested as a member of the Cyg OB2 association on the basis of 
its location within the confines of the association (Comeron et al. 
120021 1. The astrometric data on ALS 19631 and the geometry of 
the bow shock generated by this star, however, suggest that this 
runaway was instead ejected from the open cluster NGC6913 
centred ^ 3° 4 west of the star (Gvaramadze & Bomans 2008a). 

The bow shocks associated with Sk -66° 16 and Sk -68° 86 
are located, respectively, at ~ 3.5' (- 50 pc) and 2.5' (^ 35 pc) 
from the clusters KMHK268 (Fig.© and [SL63]495 (Fig.©. 
Sk-66° 16 is located in the Nil star-forming region, not far 
from our programme star Nl 1-026 (see Fig. [5}. We note the de- 
tection of a bow shock-like structure associated with one of the 
most massive stars in Nil, the ON2IIIf*: star Nl 1-031 (Evans 
et al. l20061 >. This structure is facing towards the centre of the par- 
ent association LH 10 (Fig.|7^). Interestingly, the radial velocity 
of Nl 1-031 is ^ 30 kms -1 greater than the median velocity of 
stars in Nil (Evans et al. 120061 1. which could be considered as 
indicating that this star is a runaway as well. 
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6 arcmin 



BI 253 




Fig. 4. MIPS 24/im image of the 30Doradus nebula and its surroundings. The positions of Sk-69° 206 and Sk-69° 288 and their 
bow shocks are marked by solid circles. The positions of the three programme stars, Sk-68° 137, BI253, and 30Dor016, are 
marked by crosses. The diamond point shows the position of R136. The dashed ellipsoid indicates the approximate boundary of the 
association LH 113. 



+ 

Sk-68 137 



N 



N 

+ Nl 1-026 



E 

LH 10 




Fig. 5. MIPS 24 /mi image of the Nil star-forming region with 
the programme star Nl 1-026 marked by a large cross. The 
small cross indicates the position of the ON2IIIf*: star Nll- 
031. The solid circle shows the position of Sk-66° 16 and its 
bow shock. The positions of the association LH 10 and the clus- 
ter KMHK268 are indicated by a dashed ellipse and a dashed 
circle, respectively. 



For the sake of completeness, we note also the detection 
of an arc-like nebula (Fig. [7})) attached to the candidate lumi- 
nous blu e varia ble R81 (HD 269128; Wolf et al. [1981] van 
Genderen 2001; cf. Gvaramadze, Kniazev & Fabrika 20103 and 
the 24 //m counterpart to the circumstellar nebula around the 09f 
star Sk-69° 279 (Weis et al. fl997l l. 

4. Discussion and conclusion 

The discovery of a bow shock produced by BI 237 lends strong 
support to the idea that this and other isolated massive stars in 
the field of the LMC are runaway stars (Walborn et al. 120021 



N 

E * ' 

[SL63] 495 

o 

Sk -68 86 

2 arcmin 



Fig. 6. MIPS 24yum image of the bow shock associated with 
Sk -68° 86 (marked by a solid circle). The position of the cluster 
[SL63] 495 is indicated by a dashed circle. 



Massey et al. [20051 Evans et al. [20061 l20T0l l. The young ages 
(~ 2 Myr) of BI237 and other 02-type field stars decidedly ar- 
gue for their peculiar velocities being explained by supernova 
explosions in binary systems (Blaauw 1961); the massive com- 
panion (primary) stars would simply have no time to end their 
lives in supernovae. Moreover, the high (measured or inferred) 
peculiar velocities of these stars cannot be accounted for within 
the framework of the binary-supernova scenario since it requires 
that the stellar supernova remnant (a 5 - 10 M black hole) 
receive an unrealistically high (<: 200 - 300 kms ) kick ve- 
locity at birth (Gvaramadze & Bomans 2008a; cf. Gvaramadze 
2009). The only viable alternative is that the massive stars were 
ejected in the field via dynamical three- or four-body encounters 
(Poveda et al. [19671 Leonar d & D uncan [19901 Kroupa [19981 
Pflamm-Altenburg & Kroupa 2006; Gvaramadze, Gualandris & 
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Fig. 7. MIPS 24 yum images of (a) a bow shock-like structure as- 
sociated with Nl 1-031 and (b) an arc-like nebula attached to the 
candidate luminous blue variable R81 (HD 269128). The posi- 
tions of both stars are indicated by circles. The orientation and 
the scale of the images are the same. 



Portegues Zwart 2008, 2009). Naturally, less massive (late B- 
type) stars are also ejected from their birth clusters by dynami- 
cal interactions (e.g. Kroupa 1998), but they would be difficult 
to observe in the LMC. 

The large separation of some of the 02-type field stars from 
their plausible birthplaces implies that these stars were ejected 
soon after birth (which also argues against the binary-supernova 
scenario). This implication has an important consequence for un- 
derstanding the early dynamical evolution of star clusters since it 
suggests that mass segregation in young clusters (the necessary 
condition for effective production of runaway OB stars) should 
be primordial rather than caused by the Spitzer instability. For 
example, R136 was found to already be mass-segregated at its 
age of about 2 Myr or younger (Campbell et al. 1992; Hunter et 
al. [19951 Brandl et al. [19961 de Grijs et al. 120021 . But the Spitzer 
instability could be a very fast (5> 0.5 Myr) process if the birth 
cluster is very dense (e.g. Kroupa 2008). High-precision proper 
motion measurements for the massive field stars are therefore 
required to determine the timing of their ejections, thereby dis- 
tinguishing between the primordial and the dynamical origins 
of mass segregation in young clusters (Gvaramadze & Bomans 
2008b). Future proper motion measurements with the space as- 
trometry mission Gaia will allow us to solve this problem. At 
the same time, N-body experiments are required to quantify the 
expected differences between the two types of mass segregation 
in terms of the ejection of massive stars. 

To conclude, the search for bow shocks in star-forming re- 
gions and subsequent identification of their associated stars serve 
as a useful tool for detecting runaway OB stars (e.g. Gvaramadze 
& Bomans 2008b ; Gvaramadze et al. 1201 Obi ), hence for con- 
straining the dynamical evolution of their parent clusters. Further 
search for bow shocks around young massive clusters and OB as- 
sociations in the LMC (when necessary, accompanied by follow- 
up spectroscopy of their associated stars) is therefore warranted. 
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